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ABSTRACT 

We report the discovery of a kilohertz quasi-periodic oscillation (kHz QPO) in the low-mass X-ray binary 
(LMXB) 2S 0918-549. The kHz QPO has a frequency of 1156±9 Hz, a FWHM of 70±20 Hz, and a frac- 
tional rms amplitude of 18%±2%. We also detected for the first time a type I X-ray burst in this source. The 
compact object must therefore be a neutron star. Finally, we were able to classify the source as an atoll source 
exhibiting the full complement of spectral/timing states of this class (island and banana). From the peak burst flux 
an upper limit on the persistent luminosity can be derived of 0.5% of the Eddington luminosity, making 2S 0918- 
549 one of the least luminous LMXBs showing a kHz QPO. We compare the fractional rms amplitudes of the 
upper kHz QPO across the ensemble of LMXBs. We find a strong anticorrelation with luminosity. In LMXBs 
with luminosities ^100 times lower than those of Z-sources, the fractional rms amplitude is a factor ^10 larger. 

Subject headings: accretion, accretion disks — stars: individual (2S 0918-549) — stars: neutron — X-rays: stars 



1. INTRODUCTION 

Observations with the Rossi X-ray Timing Explorer (RXTE) 
satellite have revealed the presence of several quasi-periodic 
phenomena in the Fourier power spectra of low-mass X-ray 
binaries (LMXBs). They occur on timescales similar to the dy- 
namical timescale near a neutron star, i.e. at frequencies around 
1000 Hz. Kilohertz quasi-periodic oscillations (kHz QPOs) 
have been discovered in more than 20 sources (van der Klis et 
al. 1996a,b; Strohmayer, Zhang, & Swank 1996; Strohmayeret 
al. 1996; see for the most recent review van der Klis 2000). 
■ In several LMXBs nearly coherent oscillations were discov- 
ered during type I X-ray bursts (Strohmayer et al. 1996; for 
a review see Swank 2000); these burst oscillations presumably 
occur at frequencies close to the neutron star spin frequency 
(Strohmayer et al. 1996). Recently, another high frequency 
quasi-periodic phenomenon was discovered in the power spec- 
tra of three LMXBs (Jonker, Mendez, & van der Klis 2000). 
The kHz QPOs are nearly always found in pairs, although in a 
few sources so far only one kHz QPO has been found (Zhang 
et al. 1998; Marshall & Markwardt 1999; Homan & van der 
Klis 2000). Their frequencies vary over several hundred Hz 
on timescales of hours to days; in the best studied sources, it 
is found that the frequency separation of the twin kHz QPOs, 
h>2 — v\, decreases monotonically by 50-100 Hz over the ob- 
served frequency range as the frequency of the kHz QPO in- 
creases (see e.g. van der Klis et al. 1997; Mendez et al. 1998; 
Mendez, van der Klis, & van Paradijs 1998; Mendez & van der 
Klis 1999); observations of the kHz QPO frequencies in the 
other sources are consistent with this trend (Jonker et al. 1998, 
Psaltis et al. 1998). 

Various models exist for these QPOs. Immediately after their 
discovery a beat frequency model was proposed (Strohmayer 
et al. 1996), of which the sonic-point model is the most cur- 



rent (Miller, Lamb, & Psaltis 1998; revised by Lamb & Miller 
2000). Later, Stella & Vietri (1999) proposed the relativis- 
tic precession model (but see Markovic & Lamb 2000) and 
Osherovich & Titarchuk (1999) introduced the two-oscillator 
model. In the latter model the QPO at v x (the lower kHz QPO) 
occurs at the Keplerian frequency of material orbiting the neu- 
tron star, whereas in the other two models the QPO at V2 (the 
upper kHz QPO) is the one that is Keplerian. Recently, Psaltis 
& Norman (2000) proposed the transition layer model in which 
the twin kHz QPO peaks arise due to disk oscillations at fre- 
quencies close to those predicted by the relativistic precession 
model. 

In this Letter, we report the discovery of a kHz QPO and a type 
I X-ray burst in 2S 0918-549. The source turns out to be one 
of the least luminous kHz QPO sources, however, the fractional 
rms amplitude of the QPO is quite high. We study the depen- 
dence of the fractional rms amplitudes of the upper kHz QPO 
on source luminosity across the ensemble of LMXBs and show 
that they are strongly anticorrelated. 

2. OBSERVATIONS AND ANALYSIS 

We have used observations obtained with the proportional 
counter array (PCA; Jahoda et al. 1996) onboard the RXTE 
satellite (Bradt, Rothschild, & Swank 1993). A log of the ob- 
servations and the average source count rates at the time of the 
observation can be found in Table 1 . In total ~75 ksec of data 
were used in our analysis. Data were obtained in four different 
modes. The Standard 1 and 2 modes, which are always opera- 
tional, respectively provide data with a time resolution of 0. 125 
seconds in just one energy bin and 16 seconds in 129 energy 
bins covering the effective PCA 2 to 60 keV range. Addition- 
ally, data were obtained using the mode providing the highest 
time resolution (~ 1/is), combined with the full energy res- 
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olution RXTE can provide (256 energy bins covering the 2 to 
60 keV; the GoodXenon mode). This mode saturates on count 
rates higher than —8000 counts per second. 
In order to follow spectral variations we created color-color di- 
agrams (CDs) by plotting a soft color vs. a hard color. The soft 
color was defined as the ratio between the count rates in the 
3.5-6.4 keV and 2-3.5 keV energy band; the hard color was 
defined as that between the 9.7-16.0 keV and 6.4-9.7 keV en- 
ergy band. The data were background subtracted but no dead- 
time corrections were applied (the dead-time is less than —2%). 
The colors of the 1997 observations were corrected for PC A 
gain changes, using Crab observations obtained in May 1997, 
and May 2000. Assuming the Crab is constant, observed color 
changes of the Crab between the observations of May 1997 and 
May 2000 should reflect changes of the instrument's spectral 
response. The observed changes in Crab were —6% in the soft 
color and ~1% in the hard color. The colors were converted 
to those which would have been detected by the PCA on May 
2000 (cf. di Salvo et al. 2000). Residual color errors due to the 
differences between the spectra of Crab and 2S 0918-549 are 
less than the statistical errors. 

We calculated Fast Fourier Transforms with a Nyquist fre- 
quency of 4096 Hz of data segments of 16 seconds (2-60 keV). 
The power density spectra were averaged and fitted with a func- 
tion consisting of Lorentzians (to represent peaks arising in the 
power spectrum due to QPOs), a constant to account for the 
power due to Poisson counting statistics, and an exponentially 
cutoff power law to represent the noise component apparent at 
low frequencies. Errors on the fit parameters were calculated 
using Ax 2 = 1.0 (la single parameter). The 95% confidence 
upper limits were determined using A% 2 = 2.71. 

3. RESULTS 

The color-color diagram of 2S 0918-549 (see Fig. 1; left) re- 
sembles that of an atoll source. The squares represent the aver- 
age of 256 s of data from observations 1-10, the filled circles 
are 128 s averages from observation 11, and the open circles are 
128 s averages from observations 12 and 13. Together with the 
timing properties (described below) it is clear that the squares 
represent the Island state and the open and filled circles the ba- 
nana branch of an atoll source. 

The average power spectrum of observations 1-10 combined 
(see Fig. 1, right, top) shows a strong band limited noise com- 
ponent, which could be described using an exponentially cutoff 
power law with a fractional rms amplitude of 26% ±1% (in- 
tegrated over 0.1-100 Hz, 2-60 keV), a power law index of 
0.47±0.05, and a cutoff frequeny of 21 ±6 Hz. We derived a 
95% confidence upper limit on the presence of a kHz QPO with 
a full-width at half maximum (FWHM) 100 Hz in the range of 
300-1200 Hz of -30% (rms). 

We discovered a kHz QPO at 1156±9 Hz, with a FWHM of 
68 ±23 Hz, and a fractional rms amplitude of 18%±2% (2- 
60 keV) in the average power spectrum of observation 1 1 (see 
Fig. 2). In addition, the average power spectrum of observa- 
tion 1 1 (Fig. 1, right, bottom) shows a peaked noise component 
at a frequency of 25±2 Hz, with a fractional rms amplitude of 
17%±2%, and a FWHM of 22±6 Hz. To investigate the peak 
noise component and the kHz QPO further, we also calculated 
power spectra using data in the 4.0-18.0 keV energy band. The 
values measured in the 4.0-18.0 keV energy band are the same 
within the errors, although the significance of the detection of 
the kHz QPO is somewhat higher (6.1er compared to 5.7 a in 



the 2-60 keV band). We subdivided the 4.0-18.0 keV power 
spectra of observation 1 1 in two parts. The frequency of the 
kHz QPO increased from 1 126±9 Hz in the first part of obser- 
vation 1 1 to 1218±20 Hz at the end. The source moved further 
up the banana branch, from an average hard and soft color of 
0.452±0.003 and 2.49±0.01 to 0.431±0.003 and 2.64±0.01, 
respectively. The background subtracted count rate in the 4.0- 
18.0 keV energy band, the FWHM, and the fractional rms am- 
plitude are consistent with being the same in the two parts. The 
properties of the peaked noise component were consistent with 
being the same in the two selections. 

No power spectral components were detected in the average 
power spectrum of observations 12 and 13 combined. We de- 
rived an upper limit on the presence of a 100 Hz wide kHz QPO 
of -13% rms in the range between 300—1200 Hz in both the 2— 
60 keV and the 4.0-18.0 keV energy band, and an upper limit 
of 5%-6% on the presence of a power law component at fre- 
quencies in the range of 0.1 to 100 Hz in both the 2-60 keV 
and the 4.0-18.0 keV energy band, with a power law index of 
1.0. 

During observation 11a type I X-ray burst occurred (see 
Fig. 3). This is the first detection of a burst in 2S 0918-549, and 
establishes the compact object as a neutron star. As apparent in 
Fig. 3 the burst profile was complex, with a smooth exponential 
decay with an e-folding time of 8.95±0.05 s beginning — 10 
seconds after the initial rise. Due to the very high count rates 
of more than 20000 counts per second (3 detector count rate) 
the GoodXenon mode saturated. The total count rate displayed 
in Fig. 3 is for the 0.125 s Standard 1 data which were not af- 
fected by this saturation. During the part of the decay where the 
count rate was less than 8000 counts per second spectral analy- 
sis shows cooling of the blackbody spectral component. Using 
the 16 s Standard 2 mode we derive a flux for the 16 s bin con- 
taining the burst peak of 3.5 x 10~ 8 erg cm" 2 s" 1 (2-20 keV). 
Since the mean count rate in this 16 s time bin is a factor of > 
2.5 lower than the peak count rate as measured using the 1/8 s 
bin of the Standard 1 observations, we use this factor to improve 
our estimation of the peak burst flux. We did not correct for the 
deadtime effects which become important at the count rates at 
the peak of the burst (—25%), this provides a lower limit on 
the peak flux. The pre-burst flux is less than 4.2 x 10~ 10 erg 
cm~ 2 s _1 (2-20 keV). So, the persistent emission of the source 
when it is on the lower banana branch is less than 0.5 percent 
of the Eddington luminosity. Assuming the peak burst flux to 
be isotropic and at or below the Eddington limit of 2.5 x 10 38 
erg s _1 we derive an upper limit to the distance of 4.9 kpc. 

4 DISCUSSION 

We discovered a kHz QPO at a frequency of —1 156 Hz with a 
fractional rms amplitude of nearly 20% (2-60 keV) in the X- 
ray emission of the LMXB 2S 0918-549. For the first time in 
this source a type I X-ray burst was detected, establishing the 
nature of the compact object in 2S 0918-549 as a neutron star. 
We derived an upper limit to the persistent luminosity when 
the source is on the lower banana branch of 0.005 Lsdd mak- 
ing 2S 0918-549 one of the least luminous LMXBs in which 
a kHz QPO has been observed. The spectral and timing prop- 
erties indicate that this source is a low-luminous atoll source, 
displaying the full complement of spectral and timing states (is- 
land state and a full banana branch). Similar to other kHz QPO 
sources, e.g. Mendez et al. (1999), the kHz QPO was only 
observed when the source was on the lower part of the banana 
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FIG. 1 . — Left: Color-color diagram of 2S 0918-549. The data were background subtracted, the colors of the 1997 observations were converted to the May 2000 
observations using the Crab as a reference (see text), and the burst was removed. No dead-time corrections have been applied (the dead-time is ~2%). The squares 
and circles are 256 s and 128 s averages, respectively. Typical error bars are shown for each category. The pattern resembles the Island state (squares) and the banana 
branch (open and filled circles) of atoll sources (see Hasinger & van der Klis 1989). The data in which the kHz QPO was present are represented by filled circles. 
Right, top: Rms normalized, Poisson noise subtracted average power spectrum of the observations 1-10 (squares). Strong band limited noise is present, establishing 
the identification of the Island state. The solid line is the best fit to the data with a cut-off powerlaw. Right, bottom: Rms normalized, Poisson noise subtracted 
average power spectrum of the observations 1 1-13 (circles). Banana state power law noise and a peaked noise component are present. The solid line is the best fit 
to the data with a powerlaw and a Lorentzian component. 
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FIG. 2. — Average 4.0-18.0 keV power density spectrum of observation 11, normalized after Leahy et al. (1983), showing the kHz QPO at 1156±9 Hz (6.1ct 
single trial). The solid line represents the best fit to the data. 



4 




20 40 



Time since burst start (s) 

FIG. 3. — Lightcurve showing the profile of the type I X-ray burst occurring during observation 11. The X-axis denotes the time since the start of the burst. Note 
the complex burst profile and the large increase in count rate during the burst compared with the persistent emission. The count rate is that of the three detectors that 
were operational only. 



branch. Its frequency changed from 1126±8 Hz to 1218±20 
Hz, when the source moved further on the banana branch, sim- 
ilar to other atoll sources (e.g. Mendez et al. 1999). Its fre- 
quency is similar to the frequency of the upper kHz QPO in 
other atoll sources in a similar lower-banana state (e.g. di Salvo 
et al. 2000, Mendez et al. 1999). The frequency of the lower 
kHz QPO is usually lower in that part of the CD. Furthermore, 
when both QPOs of the twin kHz QPOs are detected simulta- 
neously the frequency of the lower kHz QPO has never been 
found at frequencies larger than ^900 Hz in any LMXB (see 
van der Klis 2000). Therefore, we conclude that this kHz QPO 
is most likely the upper peak of the kHz QPO pair which is 
detected in most other LMXBs. 

After 4U 0614+109, 2S 0918-549 is another example of a very 
low luminous LMXB whose kHz QPOs are strong, whereas in 
the most luminous LMBXs, the Z-sources, the kHz QPOs tend 
to be weak (see also the review of van der Klis 2000). To further 
investigate this, we compared the fractional rms amplitudes of 
the kHz QPO in 2S 09 1 8-549 with those of the upper kHz QPO 
in other LMXBs. In order to account for the fact that in a given 
source the fractional rms amplitude of the kHz QPO changes 
as its frequency changes we measured the fractional rms am- 
plitude at both the lowest and highest frequency at which the 
upper kHz QPO is found. From extensive studies for several 
sources it is clear that the fractional rms amplitude of the up- 
per kHz QPO decreases monotonically as the frequency of the 
upper peak increases (e.g. Wijnands et al. 1997, Jonker et al. 
1998, Wijnands et al. 1998, Mendez, van der Klis, & Ford 
2000). This ensures that the method we used will yield a good 
indication of the range in fractional rms amplitude spanned in 
each source. The sample we studied here contains only those 
sources for which a significant detection of the upper kHz QPO 
in the 5-60 keV energy range was possible (4U 0614+09; Ford 
et al. 1996; 4U 1915-05, Boirin et al. 2000; 4U 1728-34, 
Strohmayer, Zhang, & Swank 1996; 4U 1608-52, Mendez et al. 
1998; 4U 1702-42, Markwardt et al. 1999; 4U 1636-53, Zhang 
et al. 1996; 4U 1820-30, Smale et al. 1996; XTE J2123-058, 



Homan et al. 1999; KS 1731-260, Wijnands & van der Klis 
1997), or in case of the Z-sources, the 5-60 keV fractional rms 
amplitude that was reported in the literature (Cyg X-2, Wij- 
nands et al. 1998; GX 17+2, Wijnands et al. 1997; GX 340+0, 
Jonker et al. 1998). Sco X-l was not included since the frac- 
tional rms amplitude of the QPOs can not be calculated accu- 
rately at present due to the severe dead-time effects present in 
this source. We calculated the 5-60 keV fractional rms ampli- 
tude in order to avoid as much as possible systematic effects on 
the derived fractional rms amplitudes of the peaks due to the 
different Njj values towards the various sources. 
The range in fractional rms amplitude of the upper kHz QPO 
in each source is plotted in Fig. 4 as a function of the simulta- 
neously measured luminosity (for the distances and the method 
used to calculate these luminosities see Ford et al. 2000 and 
references therein; the distance to 4U 1915-05 was taken from 
Smale et al. 1988). It is apparent that the upper kHz QPOs 
detected in low-luminosity sources have relatively large frac- 
tional rms amplitudes, when compared with the fractional rms 
amplitude of the kHz QPO in the bright LMXBs (e.g. the Z- 
sources). The frequencies of the end-points for each source are 
given in Table 2. 

Van Paradijs & van der Klis (1994) found that the hardness, de- 
fined as the ratio of counts between the 40-80 keV band and the 
13-25 keV band, is in neutron star LMXBs correlated with the 
X-ray luminosity. This suggests that the X-ray flux variabil- 
ity of the source at the kHz QPO frequency is correlated with 
the hardness in that band (see also Ford et al. 1997). Within a 
source a similar correlation between the fractional rms ampli- 
tude of the upper kHz QPO and the hardness of the energy spec- 
trum has been found; the harder the source the higher the frac- 
tional rms amplitude of the upper kHz QPO (e.g. 4U 0614+091, 
Ford et al. 1997, van Straaten et al. 2000; 4U 1608-52, Mendez 
et al. 1999; Aql X-l, Reig et al. 2000; 4U 1728-34, Di Salvo et 
al. 2000; Cyg X-2, Wijnands et al. 1998; GX 17+2, Wijnands 
et al. 1997; GX 340+0, Jonker et al. 1998). 
A two-component accretion flow (radially and through a disk, 
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Gosh and Lamb 1979) has been considered to explain that the 
frequency of the kHz QPO does not depend on source luminos- 
ity (Ford et al. 2000). The observed anticorrelation between the 
fractional rms amplitude of the upper kHz QPO and the source 
luminosity cannot solely be explained by an increase of an un- 
modulated part of the accretion flow towards higher source lu- 
minosity, since the increase in luminosity is too large compared 
with the decrease in fractional rms amplitude. Therefore, this 
cannot explain the findings of Ford et al. (2000). The same con- 
clusion was drawn by Mendez et al. (2000) on other grounds. 
It is important to consider the anticorrelation presented in this 



work and the absence of a correlation reported by Ford et al. 
(2000) in mechanisms producing kHz QPOs. 
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Table 1 

Log of the observations of 2S 0918-549 used in this analysis. The count rates in column 5 are the average, 

BACKGROUND SUBTRACTED, 2-60 KEV COUNT RATES OBTAINED AS IF 5 PCUS WERE ALWAYS OPERATIONAL. 





Observation 


Date & 


Amount of data 


Average count rate 




ID 


Start time (UTC) 


(ksec) 


cnts/sec (2-60 keV) 


1 


20071-11-01-02 


02-05-1997 23:57 


-4.2 


-36 


2 


20071-11-01-03 


03-05-1997 02:15 


-2.1 


-36 


3 


20071-11-01-04 


03-05-1997 06:42 


-1.8 


-45 


4 


20071-11-01-05 


03-05-1997 08:26 


-1.4 


-50 


5 


20071-11-01-06 


04-05-1997 08:27 


-1.5 


-53 


6 


20071-11-01-00 


04-05-1997 22:41 


-11.3 


-47 


7 


20071-11-01-01 


05-05-1997 22:16 


-6.0 


-48 


8 


20071-11-01-07 


09-05-1997 10:10 


-2.4 


-54 


9 


20064-06-01-00 


15-08-1997 01:38 


-4.8 


-84 


10 


20064-06-02-00 


21-09-1997 22:51 


-5.0 


-68 


11 


50060-01-01-00 


12-05-2000 12:24 


-30.0 


-135 


12 


50060-01-01-01 


14-05-2000 19:43 


-2.1 


-187 


13 


50060-01-01-02 


15-05-2000 04:59 


-2.4 


-208 
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FIG. 4. — Luminosity versus fractional rms amplitude (5-60 keV) of the upper kHz QPO in various LMXBs. Apparent is the decrease in fractional rms amplitude 
as the source luminosity increases, both within one source and between sources. 



Table 2 

The lowest and highest frequency of the upper kHz QPO for each source used to create Fig. 4. 



Source name 


Lowest frequency 


Highest frequency 




upper kHz (Hz) 


upper kHz (Hz) 


2S 0918-549 


1156±9 


1156±9 


4U 0614+09 


421±13 


1161±5 


4U 1915-05 


542±14 


1013±8 


4U 1608-52 


883±4 


1091±8 


4U 1728-34 


552±7 


1139±13 


4U 1702-42 


1000±9 


1059±11 


4U 1636-53 


1152±4 


1213±11 


4U 1820-30 


675 ±4 


1050±4 


XTE J2123-058 


1129±8 


1129±8 


KS 1731-260 


1176±3 


1176±3 


CygX-2 


731±19 


1007±15 


GX 17+2 


645 ±9 


1086±8 


GX 340+0 


535±70 


840±20 



